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We investigate the eigenvalues and eigenvectors of the staggered Dirac operator in the vicinity of the chiral 
phase transition of quenched SU(3) lattice gauge theory. We consider both the global features of the spectrum 
and the local correlations. In the chirally symmetric phase, the local correlations in the bulk of the spectrum are 
still described by random matrix theory, and we investigate the dependence of the bulk Thouless energy on the 
simulation parameters. At and above the critical point, the properties of the low- lying Dirac eigenvalues depend 
on the Z3-phase of the Polyakov loop. In the real phase, they are no longer described by chiral random matrix 
theory. We also investigate the localization properties of the Dirac eigenvectors in the different ^3-phases. 



1. INTRODUCTION 

The theoretical understanding of the Dirac 
spectrum has improved considerably in the past 
few years. Using a variety of methods such 
as finite volume partition functions H, partially 
quenched chiral perturbation theory B] , and chi- 
ral random matrix theory (RMT) || , it has been 
shown that in the phase in which chiral symme- 
try is spontaneously broken, the distribution and 
the correlations of the low-lying Dirac eigenval- 
ues are described by relatively simple universal 
functions. This description is valid in a regime in 
which the zero-momentum modes dominate the 
effective Lagrangian. The energy scale which lim- 
its this regime is known as the Thouless energy. 
For a review, we refer to Ref. 

Recently, two studies [||J|] have appeared in 
which the Dirac spectrum was investigated in 
this context for temperatures T close to the crit- 
ical temperature T c of the chiral phase transi- 
tion. This is an interesting problem, since the 
above-mentioned approach only works in the bro- 
ken phase, and one would like to find out what 
happens to the universal features as one crosses 
T c . The present study addresses these and re- 
lated questions. In addition, we also investigate 
the properties of the Dirac eigenvectors. 



* Presented by TW at Lattice 2000, Bangalore, India. 



2. Z 3 ENSEMBLES 

We are working in the quenched approximation 
using the staggered discretization of the Dirac op- 
erator. In a study of chiral symmetry restora- 
tion in the quenched approximation, an interest- 
ing observation was made in Ref. 0. The gauge 
action has a Z3 symmetry (for N c — 3 colors) 
which is broken in the deconfinement phase. As 
a result, the phases of the Polyakov loop P clus- 
ter around the elements of Z% in the complex 
plane, and one can divide the total ensemble of 
gauge field configurations into three classes with 
arg(P) = 0, ±27r/3. An example is shown in 
Fig. |l} It was found JtJ that the chiral conden- 
sate computed from the class of configurations 
with arg(P) = vanishes above T c as expected. 
However, for arg(P) = ±27r/3 it remains nonzero 
in a certain range of T above T c . This behav- 
ior can be understood qualitatively in NJL-type 
models @,| and in RMT jlO). The point is that 
the boundary conditions of the Dirac operator are 
not invariant under Z3 transformations, and for 
arg(P) = ±27r/3, the new boundary conditions 
lead to a shift of the Dirac eigenvalues to lower 
values Jio) . Using the Banks-Casher relation jll} , 
this implies a nonzero chiral condensate. 

In the following analysis, we therefore sepa- 
rate our configurations into two ensembles, those 
with arg(P) = (ensemble El) and those with 
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Figure 1. Scatter plot of the Polyakov loop in 
the complex plane in the confinement (left) and 
deconfinement (right) phase. 

arg(P) = ±27r/3 (ensemble E2), respectively. 
This separation can be done unambiguously. For 
the purpose of the present analysis, the two 
classes arg(P) = ±2ir/3 are equivalent and can 
be combined in E2. In full QCD, the fermion de- 
terminant suppresses the E2 configurations. (In 
Ref. Q, the E2-configurations were Z 3 -rotated 
before the Dirac operator was diagonalized so 
that arg(P) = in their analysis [[l2[.) 

Strictly speaking, we should distinguish three 
critical temperatures, Td for the deconfinement 
phase transition, and T c \ and T C 2 for the chi- 
ral phase transitions of El and E2, respectively. 
Here, we assume that Td ~ T c %. Since we are 
mainly interested in the region T k, T c i, we write 
T c instead of T c \ in the following. 

3. DIRAC SPECTRUM 

We have worked on 2Vf x N t lattices with Nt=4 
and 6 for which f3 c (N s -> oo)=5.6925 and 5.8941, 
respectively Jl3| . An example for the global spec- 
tral density of the Dirac operator near zero for 
T > T c is shown in Fig. ||. For the El-ensemble, 
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Figure 2. Global spectral density of the Dirac op- 
erator at T > T c for the two different ensembles. 



we find p(0) « which, according to the Banks- 
Casher relation, implies that chiral symmetry is 
restored. On the other hand, for E2 we observe 
that p(0) 7^ which implies a nonzero chiral con- 
densate, in agreement with 

We therefore expect the distribution of the low- 
lying Dirac eigenvalues to be different in the two 
ensembles. Here, we concentrate on the smallest 
positive eigenvalue, A m i n . Depending on whether 
or not (tp^) 7^ 0; t ne expectation value (A m i n ) 
scales as follows: 

T<T C : A min ~I/-\ 

T>T C : \ min ~V°, (1) 

T — T • \ • ro T/- ,5 /( <5 + 1 ) 

where 5 is one of the universal critical exponents 
of a second order phase transition. For T < T c , 
the distribution of A m i n is described by the RMT 
result, P(A mi n) = (c 2 A min /2) exp(-(cA min ) 2 /4) 
with c = V\{ipi>)\. At T = T c , RMT predictions 
for P(A m i n ) [^| are model dependent. 

Our results for P(A min ), along with fits to the 
RMT prediction in the broken phase, are shown 
in Fig. H As expected, the results for the two 
ensembles El and E2 are very different. The top 
figure corresponds to T w T c , and the bottom fig- 
ure to T slightly above T c . In both figures, (A m i n ) 
is much larger in El than in E2. This reflects the 
fact that in the symmetric phase, small eigenval- 
ues are suppressed. Also, P(A m ; n ) in El is clearly 
not described by the RMT prediction (which is 
valid for T < T c ). For E2 at T w T c , however, 
P(A m in) is still very well described by RMT, con- 
sistent with the fact that chiral symmetry is still 
broken for this ensemble. For T > T c the agree- 
ment becomes worse. 

Other quantities such as the microscopic den- 
sity could be analyzed in exactly the same way, 
and the results and conclusions will be similar. 

4. THOULESS ENERGY FOR T > T c 

The Thouless energy [ fL5f is the limiting energy 
above which the universal description of the Dirac 
spectrum in terms of RMT is no longer valid. One 
has to distinguish between the Thouless energy 
at the hard edge and in the bulk of the spec- 
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Figure 3. Distribution of the smallest Dirac eigen- 
value (from 2000 configurations) in the two en- 
sembles for T w T c (top) and T > T c (bottom). 

trum. The Thouless energy at the hard edge is 
very well understood, both theoretically |lj| and 
on the lattice ji~7| . For T > T c the hard edge is no 
longer described by RMT so that the concept of a 
Thouless energy no longer exists. However, in the 
bulk the local spectral correlations are still given 
by RMT, and it is interesting to study the bulk 
Thouless energy above T c . A convenient measure 
of the bulk spectral correlations is the number 
variance defined by 



Y?{L) = {{n{L)-{n{L))) 2 ) 



(2) 



where n(L) is the number of levels in an inter- 
val of length L after the spectrum has been un- 
folded. There are several questions related to how 
the spectrum should be unfolded, see Ref. Q for 
a comprehensive discussion. We have used en- 
semble averaging to construct the average spec- 




Figure 4. Number variance of the Dirac eigenval- 
ues in the bulk of the spectrum. Ao is the starting 
point of the interval of length L. 

tral density used in the unfolding procedure. The 
number variance was also computed by ensemble 
averaging. A typical example for E 2 (L), averaged 
over 600 independent configurations, is shown in 
Fig. |[ along with the parameter-free prediction 
of RMT. We observe that for small values of L, 
the lattice data are nicely described by RMT. We 
also see that there is a critical scale L c , the Thou- 
less scale, above which nonuniversal behavior sets 
in. In order to extract this scale from the data, 
we construct the ratio 



ratio(L) 



^rmt(^) 



^rmt(^) 



(3) 



which should start to deviate strongly from 
above L c . We find that the numerical value of 
the Thouless scale depends on where we are in 
the spectrum, i.e., on the starting point of the in- 
terval of length L. This means that spectral av- 
eraging must not be used to construct £ 2 (£) for 
the purpose of extracting the Thouless energy. 

In Fig. ||, we show the ratio of Eq. (^) for two 
different starting points, Ao = 0.5 and Ao = 0.8, 
and for various lattice sizes at /? = 6.0, which 
is above f3 c for both values of N t . Again, we 
used 600 independent configurations per parame- 
ter set. We observe that the Thouless scale seems 
to be independent of N s , but depends on N t . We 
are currently investigating the form of this Nt de- 
pendence. 
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Figure 5. The ratio denned in Eq. (^) for two 
different values of Ao and various lattice sizes. 

5. LOCALIZATION PROPERTIES OF 
DIRAC EIGENVECTORS 

In condensed matter physics, the question of 
whether or not a disordered mesoscopic sample is 
a metal or an insulator can be answered by con- 
structing the so-called inverse participation ratio 
I2, which is a measure of how many sites con- 
tribute significantly to the wave function. For the 
case of QCD, we introduce the gauge-invariant 
definition 



i 2 (x)^v ^ px{x)2 2 , 



(4) 



where V is the lattice volume and p\(x) is the 
gauge-invariant probability density 



E 

a=l 



(5) 



Here, x denotes a lattice site, a is a color in- 
dex, and ipx( x ) 1S a component of the eigenvector 
corresponding to eigenvalue A. Because of chiral 
symmetry we have J2 X cvcn P\(x) = Ex odd Pa (a;)- 
A completely localized state therefore has I2 = 
V/2. This case corresponds to uncorrelated eigen- 
values and is described by the Poisson ensemble. 
(The corresponding mesoscopic sample would be 
an insulator.) For a completely delocalized state, 
P\{x) is the same for all x, and I2 — 1. In RMT, 
we find 



N C V 
N C V + 



1 



h (6) 



This case corresponds to extended wave functions 
with many significant components and strongly 
correlated eigenvalues. (The corresponding meso- 
scopic sample would be a metal.) In QCD, a sim- 
ple argument (see 0) shows that if chiral sym- 
metry is broken, the Dirac wave functions must 
be extended. 

Our results for the inverse participation ratio 
of the low-lying Dirac eigenvectors are shown in 
Fig. [| Several features are worth noting: 

1. Below T c all eigenvectors are extended, which 
is consistent with the fact that the eigenvalues 
are described by RMT both at the hard edge 
and in the bulk of the spectrum. The eigen- 
vectors corresponding to smaller eigenvalues are 
slightly more "localized" than those correspond- 
ing to larger eigenvalues. The data agree well 
with the RMT prediction, (I2) — 4/3, for large 
(A) . In contrast to the observation of Ref . JTjJ for 
Wilson fermions, we did not find signs of strong 
localization. 

2. The eigenvectors corresponding to eigenvalues 
in the bulk of the spectrum remain extended for 
El and E2 at all temperatures considered, con- 
sistent with the fact that the eigenvalues in the 
bulk continue to be described by RMT above T c , 

3. Most interestingly, the eigenvectors of the El 
ensemble corresponding to the lowest eigenvalues 
become more localized at and above T c , while in 
the E2 sector the lowest eigenmodes remain ex- 
tended. 

We are currently studying the topological prop- 
erties of the low-lying eigenvectors and will pre- 
sent our results in the near future. 
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Figure 6. Average inverse participation ratio (I2) 
of the low-lying Dirac eigenvectors vs average 
eigenvalue in the different ^-phases for temper- 
atures below, at, and above T c . 

6. CONCLUSIONS 

The behavior of the Dirac spectrum at and 
above the critical temperature can help to im- 
prove our understanding of the chiral phase tran- 
sition. In particular, the properties of the Dirac 
eigenvectors deserve further study. Of course, the 
quenched approximation and the use of staggered 
fcrmions make it difficult to establish contact 
with continuum QCD. These systematic prob- 
lems can be circumvented by using dynamical 
Ginsparg- Wilson fermions which, however, re- 
quire large computational resources. For the mo- 



ment, we are restricted to exploratory studies 
such as the present one. 
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